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SUMMARY 

vj'he  problem  of  acoustic  radiation  of  panels  excited  by  random  pressure 
flucluration  of  the  turbulent  boundary  layer  further  investigated.  In 

the  past.  oxperi mental  investigators  have  eoneenlrated  only  on  certain  phases 
of  the  general  problem:  for  example,  the  statistical  behavior  of  the  pressure 
field  (Hof.  1),  or  the  motion  of  the  panel  (Kef.  2<»).  Recently  the  author  presented 
a  comprehensive  set  of  experimental  results  (Kefs,  d  and  f>)  that  include  all 
statistical  information  concerning  both  the  forcing  function;  that  is,  pressure 
fluctuation  and  the  response  functions  of  the  panel  motion  as  well  as  the  radiation 
field. 

The  main  purpose  of  this  paper  is  lo  show  by  using  a  relatively  simple  functional 
representation  of  the  space-lime  correlation  of  the  wall  pressure  fluctuation, 
and  by  the  use  of  Lyons -Dyer  method  (KelV-20)^  that  motion  and  radiation  intensi¬ 
ty  of  a  simply-supported  panel  agree  reasonably  well  with  ihe-atiihor's 
experimental  results. 

The  most  striking  feature  of  the  excitation  mechanism  is  the  so-called  eoin- 
eideaee  which  has  profound  effects  on  the  response  of  the  structure  and  power 
radiations.  If,  under  certain  conditions,  a  mismatch  occurs  between  wave  speeds 
on  the  panel  and  the  pressure  field,  panel  displacement  and  acoustic  radiation 
should  be  reduced.  Such  a  mismatch  is  caused  by  a  turbulence  pressure  eddy 
which  decays  faster  than  the  mode's  wavelength  on  the  structure. 
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NOMENCLATURE 

-  plate  model  damping 

a .  b 

-  lengths  of  panel  sides 

A.  K.  N 

-  constant 

1) 

-  heading  stiffness 

E 

-  Young's  modulus 

*  f  or  w 

-  frequency  (to  -  2rrf) 

h 

“  panel  thickness 

K 

-  wave  number 

m,  n 

-  mode  number 

M 

-  Mach  number 

Jft 

-  plate  m'ass 

/T 

v  p 

-  rms  wall  pressure  fluctuations 

p(  ) 

power  spectral  density  of  the  wall  pressure  fluctuations 

PWl, 

-  total  sound  power  level 

m  ) 

-  correlation  coefficient 

ii 

-  free  stream  velocity 

l,c 

-  convection  velocity 

x.  y 

-  coordinates 

Y 

-  panel  displacement 

e 

-  eddy  lifetime 

0 

-  wave  form 

6 * 

"  boundary  layer  displacement  thickness 

r  1  -  |' 

-  time  delay 

7]  y  -  y' 

-  lateral  partial  separation 
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4  y  -  y'  “  Io«|* lluclinsil  partial  separation 
r  -  wall  shear  stress 

a  -  standard  deviation 

ft  r,  ft  ,  -  acoustic  and  structural  damping 

AC-/  f  b  1 


r, 


m ,  n 


in .  n 


-  eigenvalue 

-  total  damping  ratio 


SHEEI  10 


AC 


0-7O0( 


r 


NUMBER D6- 99^  “  VGL.III 
REV  LIR 


TKST  RKSUl.TS  FROM  IU tUNDAKY  LAVKIt  FACILITY 
(THKOKY  AND  KXPFKIM  KNTAl.  COM  PARISON) 

I.  INTRODUCTION 

The  exterior  surface  of  an  aircraft  is  excited  directly  by  turbulence  as  well  as  by 
sound  generated  by  the  turbulent  boundary  layer.  The  aircraft  skin  is  flexible  and 
the  fluctuating  pressure  in  the  boundary  lajer  (forcing  field)  causes  the  surface  to 
vibrate.  This  surface  vibration  acts  as  a  radiator  of  sound.  Measurements  of 
the  wall  pressure  flucturntion,  panel  i-espon.se,  and  acoustic  radiation  character¬ 
istics  of  typical  panels  have  been  reported  by  many  investigators.  The  intention 
in  this  paper  is  to  discuss  some  of  these  measurements  (primarily  those  made  by 

the  author)  and  to  develop  a  theoretical  model  for  comparison  with  experiments. 

Consider  the  surface  pressure  fluctuation,  since  the  fiuclurution  force  on  the 
wall  in  any  frequency  band  causes  the  panel  structure  to  vibrate  at  wavelengths 
of  the  order  of  the  scale  of  the  pressure  field.  Typical  wall  pressure  power 
spectra  and  correlation  fields  can  be  seen  in  Itefs.  1  through  7  for  subsonic  flow 
and  in  Kef.  H  for  supersonic  flow.  These  data  show  that  the  mean  square 
pressure  fiuclurution  is  a  slowly  varying  function  of  Mach  number.  These 
spectra  are  in  good  agreement  for  Struhal  numbers  u>6f/U  >0.2.  At  lower 
values,  the  various  measurenjcnls  do  not  agree  well.  The  differences  in  results, 
may  possibly  be  due  to  variations  in  Reynolds  number  or  Mach  number  or  to 
extraneous  noise.  It  is  not  possible  to  find  the  characteristic  length  and  velocity 
that  eliminates  the  Reynold's  number  effect.  Reference  !)  shows  for  instance 

rr> 

that  the  ratio  ^  p  /T  approaches  a  constant  only  lor  very  high  Reynold's 
w  t, 

number.  Furthermore,  the  non  dimensional  power  spectrum  p(w)U/r  “g* 

shows  similar  behavior  only  in  a  limited  nondimensional  wave  number  range. 

The  Mach  number  range  of  the  various  experiments  extend  from  zero  to  almost 

five,  and  it  is  not  surprising  that  differences  in  the  shape  of  the  spectrum  are 
found. 
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Interference  al  low  frequencies  due  to  extraneous  noise  sources  is  commonly 
encountered  in  wind  tunnel  facilities.  The  cause  of  this  noise  in  the  various  types 
of  wind  tunnels  depends  on  type  of  blower  system  and  on  the  air  control,  and 
installation  environment.  In  the  present  test  facility.  1 5  feel  of  acoustic  muffler 
are  used  ahead  of  the  settling  chamber,  Kaeh  component  is  mounted  on  an 
independent  concrete  pad  floating  on  sand.  However,  interference  still  occurs  at 
frequencies  of  100  cps  and  lower  and  varies  somewhat  with  the  flow  speed  and 
diffuser  setting  (Ref.  -1).  K  vidence  of  acoustic  interference  is  shown  by  the 
narrow  band  spatial  correlation  measurements  of  the  wall  pressure  made  at  low 
frequencies  (Ref.  10).  These  measurements  clearly  show  a  very  strong  acoustic 
pressure  superimposed  on  the  turbulent  pressure.  It  may  only  be  practical  to 
limit  low  frequency  power  spectrum  measurements  to  flight  tests. 


The  response  of  panels  to  turbulence  lias  also  been  the  subject  of  various 
papers  from  both  experimental  and  theoretical  aspects  (Refs.  11  to  20).  These 
papers  are  concerned  primarily  with  the  vibration  correlation  response  near 
the  free  bending  wave  length  of  the  panel,  which  accounts  for  most  of  the  excita¬ 
tion.  Theory  in  Ref.  Hi  is  based  on  the  acceptance  of  the  pressure  field  by  the 
beam,  indicating  the  influence  on  the  displacement  response  and,  hence,  the 
stress  of  the  structure.  These  /  esenaiU  effects  of  matching  the  bending  wave 
number  and  frequency  of  the  panel  with  the  turbulence  wave  number  and 
frequency  (aero-dynamic  coincidence)  have  been  observed  experimentally 
(Refs.  f>,  20,  21,  and  22).  Reference  f>  extends  (he  experimental  results  lor 
the  conditions  below  and  above  coincidence  both  in  terms  of  model  displacement, 
panel  stresses  and  acoustic  radiation.  Reference  21  includes  a  brief  resume  of 
the  behavior  of  the  panel  displacement  obtained  by  space-time  correlation  over 
a  complex  panel.  These  results  are  associated  with  the  measured  total  mean 
square  displacement  and  total  acoustic  power  radiated. 
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It  turns  out  that  the  amount  of  energy  fed  into  the  structure  will  depend  on  the 
wave  number  component  and  amplitude  of  the  pressure  field.  On  the  other  hand, 
the  amount  of  energy  accepted  by  the  structure  will  depend  on  how  many  struc¬ 
tural  modes  the  pressure  lield  will  excite.  However,  not  all  the  energy  accepted 
by  the  structure  will  be  converted  into  displacement  amplitude  but  will  be 
partially  dissipated  internally  and  partially  radiated  back.  This  loss  corre¬ 
sponds  to  the  total  damping  consisting  of  acoustic  and  sfructural  damping. 

Finally,  the  acoustic  power  radiated  has  been  tlie  subject  of  considerable 
investigation,  expecially  by  the  Bolt,  Beranek,  and  Newman  Corporation 
(Kefs,  25  to  27).  Their  analysis  has  shown  that  the  sound  power  radiated  from 
a  reverberant  field  for  simply-supponed  panels  is  proportional  to  the  perim¬ 
eter  of  the  panel  for  vibration  up  to  the  critical  frequency.  When  comparing 
simply-supported  or  (damped  boundaries,  it  was  found  that  the  clamped  case  was 
the  most  efficient  radiator.  Measurements  of  (he  acoustic  damping  for  various 
modes  of  vibration  have  been  made  and  reported  by  the  author,  (Kef.  -1).  The 
change  in  radiated  power  with  velocity  has  shown  the  effect  of  aerodynamic 

coincidence,  as  well  as  a  marked  cancellation  effect  which  reduces  the  acoustic 
5  2 ,  .'I 

power  law  from  a  M  to  a  M  '  dependence.  The  change  is  caused  by  a  pro¬ 
gressive  increase  in  wave  number  for  coincidence  where  a  similar  wave 
number  of  the  panel  mode  becomes  larger  than  the  corresponding  wave 
number  of  the  turbulent  pressure  field. 

Owing  to  the  complexity  of  the  problem,  a  measured  input  covariance  was  used 
in  predicting  panel  response  due  to  turbulence.  The  panel  is  flat,  rectangular 
and  has  either  simply-supported  or  clamped  edges.  The  vibration  of  the 
surface  is  coupled  to  the  turbulent  pressure  field  beneath  the  panel  and  to  the 
radiation  field. 
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II.  PROPERTIES  OF  THE  WALL  PRESSURE  FLUCTUATION 
We  have  proposed  in  the  previous  work  Ref  (5)  two  functional  representations  of 
the  w all  pressure  covariances  based  on  experimental  observations  of  data  from 
a  turbulent  channel  flow.  Since  the  structural  panel  is  excited  randomly,  i.e.  by 
turbulence,  the  How  model  is  useful  in  studying  coupling  between  the  wall 
pressure  and  the  panel  surface.  The  model  consists  of  a  con  reeled  wall  pros 
sure  pattern  with  space  and  time  dependence.  It  retains  the  characteristics  of 
the  superposition  of  a  wave  system  with  phase  and  amplitude  related  to  wave 
number  and  frequency  spectrum.  This  model  may  not  be  sufficiently  rigorous 
to  insure  a  detailed  description  of  the  flow  field,  but  it  proves  to  be  us  'ul  m 
describing  the  response  of  the  panel  structure  to  turbulence  because  it  $f» ovules 
a  means  for  studying  the  effect  of  pressure  fluctuations  on  the  motion  of  the  panel. 


The  representation  chosen  here  to  describe  the  wall  pressure  tlueitialion  in  a 
moving  frame  is  given  by  the  linear  combination  of  two  Gaussian  distributions .  each 
of  which  has  a  maximum  value  that  decays  in  time  and  changes  in  separation  with 
time.  The  behavior  in  the  moving  frame  reference  of  the  cross 'Correlation 
function  U(£,i}.t)  has  been  represented  by: 


<lXx,y.l)p(x-i-g,yi  ij.ti t)> 


ITT  -lu-VWI/a,,2 
1  V  'A., 


After  decomposing  each  individual  correlation  H(f  .r)  into  (wo  Gaussian  distribu¬ 
tions.  an  "upper"  one  and  u  "lower"  one.  a  single  decay  rather  than  an  individual 
decay  for  each  of  IhcGuussinn  distributions  w  as  used.  The  doca\  of  the  cross - 
correlation  is  given  by  exp  (|£|  /Uc0)*  rather  than  by  e.\p  (|r|  0).  as  in  the  previous 


'Appendix  A  discussed  liu*  clmne  using  one  ul  l\vu  tonus  ol  decay 
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paper.  The  panel  does  not  respond  similarly  for  all  frequencies  and  wave  numbers 
uf  the  turbulence;  certain  combinations  of  wave  number  and  frequency  are  excited 
more  strongly  than  others.  When  the  wall  pressure  matches  both  the  wave* 
number  and  tiie  frequency  of  a  particular  panel  mode,  a  condition  called  coincidence 
occurs.  This  results  in  a  very  strong  excitation  of  running  flexural  waves  on 
thcpancl.  Then  the  amplitudes  of  any  Fourier  components  of  the  panel  displace¬ 
ment  depend  on  the  amplitudes  of  the  characteristic  lengths  of  the  pressure 
fields.  U  $ .  The  characteristic  length  or  correlation  length  is  a  function  of 
both  wave  number  and  frequency,  implying  that  U  (Kj  ,u>).  Tlu*  moving  frame 
spectrum  shows  that  a  very  uniform  energy  content  up  to  the  wave  number 
Kj  I  /It  0  exists,  beyond  which  the  energy  diminishes  rapidly.  The  following 
spectrum  was  obtained  from  a  Fourier  transform  of  the  envelope  of  the  correlation 
maxima  exp  -Ul/U  ,0.  Figure  1  shows  the  eddy  lifetime  and  the  moving  frame 

spectrum. 


1><W 


U  0 

HH1  !Ki  ucr> 


'I'he  e<|iii  valent  representation  in  terms  of  the  wave  number  and  frequency  spectrum  is 
obtained  from  the  triple  Fourier  transform. 


I  m  -i  K  «IK.,1PWT) 

1>(K  K.;,w) - --JjJm.’UT)  e  1  "  d^dydT 

(2tt)'  -oc 
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The  problem  is  to  integrate  (.*{) .  I  .el  us  first  obtain  Ironi  (I)  the  power  spectrum. 


Vl 

P(0,0,w)  —  e 


•>  •>  •>  •}  9  •> 

a.,,.,  _V“‘/2lv 


u(.  />» 


11  y-2» 

e  v 


ami  the  autocorrelation 


(t,uT2/2aJ)“ 


»(0,0,r) 


A  |  e 


A2  * 


(I) 


(•r») 


The  power  spectrum  of  equal  ion  (I)  is  plotted  in  Kin.  2  and  compared  with  the-  measured 
spectrum.  A  good  lit  is  obtained  at  the  lower  frequencies  ior  u6*/U  <  1.5; 
at  the  lusher  1  requeue ics  the  Gaussian  spectrum  distribution  lulls  ol  much  luster 
than  the  measured  spectrum.  These  dillerenees  in  the  power  spectrum  were 
discussed  previously  in  ltd.  5  where  an  exponential  space-tune  correlation  was 

obtained  containing  higher  I  requeue}  components  than  the  model  described  b>  l-.q.  I  . 
Also  a  comparison  is  made  in  Kig.2  between  the  measured  autocorrelation  and 
the  one  described  111  equation  (5). 


Kroin  equations  (I)  and  (.'!).  the  longitudinal  cross  spectra!  dens  it  \  becomes: 


l*(w) 


w/l  1  1  K 

e  I 


l/V  0 
e 


]  /n  0 
<• 


(I  /ti(fl)“i(u>  i!^>K1)“  (J. 'V^0)~i(u  T'^-Kj)* 


r,  ' 


uj/H  Is, 
e  ! 


(I /II  0)“  »  {w/Ut.  '  K,)“  (I  Ml  .0)*  1  (ta'Ut.  K  f  )' 


i/0 


0i> 
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Since  l lie  cross  spectral  density  is  complex,  giving  information  about  the 
amplitude  and  phase,  equation  (<i)  can  also  be  written  as 


P(Kj  0,w)  -  l>(K, .  0.w)|e 


~*0( K j  w) 


Then  the  phase  angle  is  given  by 


0(K.  0,w)  tan  P/ot 

*  ♦ 


Similarly,  the  lateral  spectral  density 


Act  2  -£rl2a,2/2U  2  ~a\  K2*/2 

J,((VVw>  -  liir-0  c  e 

e 


Act  "°2W  /2Ue  ~a2  K2  /2 

1 1HT  0 

C 


On  combining  (ti)  and  (8).  one  obtains  the  two  dimensional  cross  power  spectrum 
P(Kj,  K„,  w)  between  the  vectors  K(  and  K.,. 

The  cross  power  spectral  density  shows  the  rale  of  transfer  of  energy  from 
small  to  large  wave  numbers.  The  most  coherent  pattern  of  the  pressure  eddy 
moves  with  the  convection  velocity  at  different  frequencies  and  wave  numbers. 

G 

Useful  definition  of  the  convection  velocity  has  been  given  Ir,  Kefs.  12  and  28, 
where  it  is  show'll  that  the  integral  time  scale  is  the  maximum  in  the  frame  of 
reference  moving  downstream  with  this  velocity.  Since  only  the  lateral  correla¬ 
tion  shows  no  moving  axis,  exp-(»l/U  0)  -  1;  the  time  scale  0  is  constant  and 
very  small  along  ij  axes.  In  all  other  components,  the  value  of  0  is  large. 

The  eddy  lifetime  0  were  assumed  to  correspond  to  the  time  in  which  the  value 
of  the  correlation  coefficient  obtained  from  the  envelope  of  the  correlation 
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maxima  dropped  to  l/e.  Since  it  is  known  that  the  convection  velocity  is  not  the 
same  for  different  wavelengths  and  frequencies,  it  is  convenient  to  define  a 
convection  velocity  o>0/K*  =  Kj)sueh  that  when  Kj  -  K*,  then  |P{Kj  0.  0)|“ 


is  a  maximum. 


The  magnitude  of  the  spectrum  given  by  Kq.  (>  is: 


2  2  2  •> 

I  U  0  (i  1  u  0  ) 


»h2 -i£*bs3  — 

'  (2 7r) "*  1  IOJ~0“(J  iUeKJ/co)“  JifcT0“(|-UcK 


By  maximizing  j  P(K.  ,  0,  w)  |  “  with  respect  to  K  where  U  _  and  0  are  constant. 
1  10 

We  observed  the  following: 

C  |P(w)|  2(‘l Uc“0“);  at  K*  -  0  when  caO-M) 


|  P(w)  |  “ - 5 - 5 ;  at  K  4  --  0  when  u>0<i 

(i/u cor  *  (w/ucr 


|P(K.  ,0,  ta)  | 

1  '  1’  ’  7 1  max 


0<K j <  00 


•->  *) 
21J  “<T 


|P(cu)  |  “  — ;  at  K*  -  0  when  c oO 
(■>*)- 


ii’(w)i2  u  .V 

uo  vLV2- 1 


1  “ Vc 0  /II  -1/U  0  <  ,r  •  and  when  o>0>l 

c  e  u 

c 


From  (10)  the  convection  velocity  UJ  is  given  by: 


IH  U  ,  / 

l'  l  V0V- ! 
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The  first  three  conditions  simply  mean  that  the  notion  of  convection  velocity  is  not 
a  useful  concept.  It  is  of  considerable  significance  that  in  the  region  w0»l  the 
peak  amplitude  varies  like  Kj1  sw/U^.  in  this  ease  UJ.  It  can  he  seen 
that  the  convection  velocity  increases  rapidly  with  decreasing  wavenumber  until 
its  speed  relative  to  the  free  stream  velocity  will  be  larger  than  the  sound 
velocity.  This  region  corresponds  to  the  Mach  wave  radiation  discussed  in  Refs. 

29  through  92. 

A  plot  of  the  measured  (real  part)  one-dimensional  wave  number  spectrum  is 
shown  in  Fig.  3.  it  is  compared  with  Uie  real  part  of  Kquation  (li).  The  mea¬ 
sured  real  part  is  obtained  by  the  Fourier  transform  of  the  narrow  band  spatial 
correlation  which  also  gives  the  phase  information.  The  differences  in  the 
lower  frequency  region  are  due  to  variation  in  U*;  at  high  frequencies  the  dif¬ 
ferences  are  introduced  by  a  measuring  technique  which  uses  a  wider  effective 
bandwidth  filter. 

Figure  4  shows  the  computation  of  the  real,  imaginary,  and  absolute  value  of 

the  longitudinal  spectrum.  The  phase  information  is  obtained  from  the  ratio  of 

the  real  and  imaginary  part.  This  figure  shows  that  energy  decreases  with 

» 

increasing  frequency  and  wave  number.  Furthermore,  it  can  be  seen  that  in 
the  lowest  wave  number  region  there  exists  a  noticeable  difference  between  the 
maximum  of  the  absolute  value  of  the  spectrum  and  that  of  the  real  part. 

Figure  5  clearly  shows  that  the  average  convection  velocity  corresponds  to  that 
range  of  wave  numbers  where  the  disturbances  are  coherent  (U  *  -  constant 
independent  of  wave  number).  On  the  other  hand,  the  lowest  wave  number  distur¬ 
bances  are  incoherent,  moving  with  significantly  different  convection  velocity. 
The  results  summarized  in  Fig.  5  confirm  the  conjecture  of  Ref.  32  that  the  low¬ 
est  wave  number  disturbance  can  produce  Much  wave  radiation  even  in  a 
subsonic  boundary  layer.  Figure  6  shows  the  region  of  supersonic  and  subsonic 
wave  velocity. 
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III.  MEASUREMENT  OK  THE  PH  ESSUH E  FLUCTUATION 


Either  static  (ires sure  fluctuation  at  a  wall  or  total  pressure  in  the  stream 
have  been  measured.  However,  measurement  of  the  stream  total  pressure 
fluctuation  in  a  channel  flow  has  been  associated  with  the  measurement  of  the 
wall  pressure  fluctuation.  The  total  pressure  fluctuation  is  useful  in  studying  the 
impingement  effect  of  the  turbulent  flow  on  the  surface.  It  is  also  presently  used 
to  measure  the  radiated  acoustic  energy  from  the  boundary  layer  by  mounting  the 
transducer  in  the  unperturbed  region  ol  the  free  stream. 

The  sensing  element  in  the  total  pressure  transducer  is  a  solid  piezoelectric 
ceramic  cylinder  with  an  effective  diameter  of  0.0-1  in.  It  is  mounted  on  a  cone 
tip  of  a  :l/l fi-inch  diameter  lube  (Fig,  7).  The  signal  output  due  to  vibration  of 
the  transducer  when  mounted  on  the  duct  was  lound  to  lie  very  small  compared  to 
the  pressure  signal.  The  vibration  increase  us  the  transducer  was  moved  into  the 
stream;  however,  up  to  an  extension  of  1 .5  inches  from  the  wall,  the  vibration 
output  was  considerably  lower  than  the  signal  due  to  pressure  fluctuations. 

To  facilitate  the  comparison  of  the  wall  pressure  with  the  stream  pressure 
fluctuations,  the  stream  pressure  transducers  were  made  to  be  adaptable  for 
flush  mounting  in  the  duel  wall.  It  was  acertamed  that  the  measured  wall 
pressure  corresponds  to  the  power  spectrum  distribution  measured  previously 
with  different  transducers.  As  expected,  the  root  mean  square  stream  pressure 
fluctuation  is  higher  than  the  wall  pressure  (Fig.  8) .  The  distribution  of  the 
pressure  across  the  duct  shows  a  maximum  near  the  wall  m  the  region  of  higher 
shear  fluctuation.  The  pressure  level  decreased  towards  the  center  ol  the  duel, 
indicating  a  behavior  similiar  to  that  of  the  velocity  fluctuation.  Figure  It  is  a 
plot  of  the  stream  pressure  spectrum  and  wall  spectrum.  The  stream  pressure 
spectrum  has  a  shape  similar  to  the  usual  velocity  spectrum. 
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Space-lime  correlation  measurements  of  the  wall  pressure  fluctuation  were  made 
over  the  fuselage  skin  of  the  727  airplane  to  determine  the  wall  pressure 
covariance  to  be  used  for  structural  and  acoustic  response  study.  The  power 
spectrum  distribution  with  frequency  is  similar  to  the  one  obtained  previously, 
Reference  -1,  on  the  7U7  prototype  airplane.  Considerable  attention  was  given 
in  the  measurement  of  the  lowest  frequency  region  in  an  attempt  to  clarify  the 
argument  about  the  behavior  of  the  power  spectrum  in  this  range. 

Rive  Bruel  and  Kjaer  Type  II.'IG,  0.25-inch  diameter  micophoncs  were 
mounted  on  !/•!— inch  thick  plates  at  various  distances  apart  in  the  direction  of 
flow.  The  plate  was  flush  mounted  on  a  cutout  portion  of  the  airplane  panel 
above  the  wing  level  approximately  71  feel  from  the  nose.  In  the  vicinity  of 
the  transducers,  a  Preston  tube  was  mounted  on  the  plate  for  the  local  skin 
friction  measurement. 


The  wall  pressure  spectra  was  recorded  for  frequencies  as  low  as  20  ops  and 


up  to  20,000  eps.  A  typical  spectrum  is  shown  in  Fig.  10  at  Mach  number  0.77. 
This  power  spectrum  does  not  rolloff  with  decreasing  frequency  like  the  Hodgson 
spectrum  measurement  over  the  wing  of  a  glider  (Ref.  2)  for  very  low  speed. 

The  differences  between  spectrum  on  the  727  airplane  and  the  one  obtained  in  the 
Boeing  boundary  layer  facility  occurs  at  a  lower  frequency  which  shows  that  the 
duct  spectrum  slightly  increases  in  amplitude  while  the  airplane  spectrum 
remains  constant. 


Figure  11  shows  that  the  ratio  of  the  root  mean  square  pressure  fluctuation  with 
the  wall  shear  stress  plotted  in  terms  of  Mach  number  for  the  727  airplane 
compares  with  the  previous  measurement  in  the  duel  flow  and  with  the  lioundary 
layer  flow  of  Refs.  1  and  7.  As  expected,  the  727  data  are  slightly  lower  than 
the  wind  tunnel  spectra  because  of  bamlwilh  limitation  of  the  instrumentation. 
However,  the  overall  picture  shows  the  trend  of  the  variation  of  the  ratio 
J  jl/Vw  with  Mach  number  lor  relatively  high  Reynold's- number. 
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Since  the  study  of  the  structural  response  to  turbulence  requires  not  only  the 
mean  square  value  of  the  pressure  field  but  also  the  energy  contained  in  a 
specific  frequency  band  corresponding  to  a  resonant  frequency  of  the  structural 
panel,  the  filter  center  frequency  and  the  bandwidth  has  to  be  selected  to  fit 
the  particular  mode  of  interest.  A  narrow  band  special  correlation  muusure- 
ment  of  the  wall  pressure  reported  in  Hef.  5  has  shown  that  the  eorrclation 
coefficient  has  a  unique  dependondoney  on  u>£/Ut  and  it  decays  as  oxp-(|£|/U.0) 
for  the  longitudinal  correlation  equation  obtained  from  (1)  after  taking  the  trans¬ 
form  with  respect  tor, 

I 

the  longitudinal  cross  power  spectrum  becomes: 


I>  (1,0,0) 


U  Of  . 
l’(w)  e  0  cos-5 
c 


-  I  Sill  -j-p 

'  c 


The  real  part  P(4,o,u){l  is  identical  to  experimental  data  obtained  by  a  narrow 
band  spacial  correlation.  The  imaginary  part  can  be  obtained  experimentally  by 
shifting  the  phase  of  one  signal  90-degrees  at  frequency w  before  correlating 
the  two  signals. 

Assuming  that  the  tiller  characteristics  are  constant  within  (w  j  -  Aw  to 

Wj  1  Aw)  and  zero  elsewhere,  the  amount  of  energy  per  bandwidth  is  related  to 

the  amount  of  energy  per  cycle  by: 


,w  t-Aw 

j  P(|,o,w)dw  HAu)  (|,o) 
wpAw 


where 


P(4,o,w)  --  7j J J_  K(4,o,T)e"IWTdT 


The  measured  spectrum  obtained  Innn  equation  (lit)  can  be  compared  with  the 
single  frequency  spectrum  of  either  equation  I:i  or  H. 
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IV.  THE  STRUCTURAL  RESPONSE  TO  A  TURBULENT  BOUNDARY  LAYER 


Since  the  pressure  field  over  the*  panel  is  responsible  for  generating  the  motion 
of  the  panel,  which  consists  of  a  multimode  vibration,  a  Fourier  analysis 
in  two  space  dimensions  £  and  r]  and  one  time  t  dimension  would  determine  how 
extensive  is  the  influence  of  the  pressure  field  in  controlling  the  response  of  the 
panel.  This  response  would  not  bo  uniform  since  certain  combinations  of  wave 
numbers  and  frequencies  in  the  component  of  (he  surface  displacement  are 
excited  more  than  others  (Ref.  5).  The  phase  relationship  for  a  particular 
frequency  and  wave  number  remains  fixed  so  that  the  resultant  excitation  will 


continually  grow  causing  *a  buildup  of  displacement  wave  amplitude  inversely  to 
the  panel  damping.  The  panel  lends  to  respond  with  runningwave  ripples  having  trace 
speeds  equal  to  the  eckly  convection  speeds  of  the  turbulent  boundary  layer 

(Refs,  ill  and  M). 


a.  Response  of  Simple  Panel 

The  vibration  response  of  simple  panels  excited  by  turbulent  boundary  layers 
is  considered,  ’rite  panel  loading  is  described  by  one  of  two  models  of  the 
pressure  covariances  obtained  from  experimental  data  from  turbulent  channel 
flow  (Ref.  5).  Numerical  computation  of  the  response  for  simply  supported 
panels  has  been  made  (Ref.  2J)  using  an  idealized  model  (Refs.  J5  and  22)  with 
pressure  covariances  having  a  delta  function  frequency  spectrum  with  exponen¬ 
tial  decay.  Since  the  results  are  qualitatively  descriptive  of  the  panel  behavior 
in  a  convcclcd  field,  it  was  fell  that  the  improvement  of  this  model  by  an 
accurate  assessment  of  the  pressure  field  may  have  a  numerical  importance 
on  the  prediction  of  the  structures  response  to  turbulence. 


Since  relatively  large  numbers  of  modes  are  present  in  a  simple  panel,  the 
response  behavior  in  each  one  of  the  modes  depend  on  the  flow  characteristics 
(u>,  Iv,  Ut>,  and  0)  and  on  the  panel  characteristics  (a,b,h,aHj  and  E).  For 
certain  combinations,  only  a  few  modes  dominate  the  sped  nun;  for  others,  the 
overall  response  is  controlled  by  relatively  large  numbers  of  modes.  A 
practical  example  is  given  by  7  x  12  x  .  01  inch  panels  excited  by  flow  Mach 
number  M  -  0.52  with  the  boundary  layer  displacement  thickness  6*  -  0.  155". 
Results  show  that  above  the  Glh  mode  (m  -  (j,  n  -  1)  the  response  is  so  small 
that  any  contribution  to  the  total  response  is  negligible.  Instead,  for  a  longer 
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panel  :1G  x  G-l/2  x  0.  01  inches  at  the  same  flow  environment,  the  mean  square 
of  the  20th  mode  (m  =  20,  n  -  1)  is  only  2  decays  below  the  first  few  modes. 

It  is  difficult  and  extremely  lengthy  in  such  situations  to  compute  the  response 
of  each  single  mode.  A  model  that  is  applicable  for  a  large  number  of  degrees 
of  freedom  becomes  very  useful. 

1 

For  a  practical  airplane  structure,  a  panel  is  divided  into  bays  of  panels  with 
a  frame  stiffener  around  the  edges  and  stringers  along  opposite  edges. 
Preliminary  measurements  show  that  the  response  of  such  structures  behaves 
as  a  running  wave  ripple,  which  favors  the  type  of  deformation  corresponding 
to  higher  order  modes.  Each  adjacent  panel  becomes  independent.  Previous 
work  on  multiple  panels  given  in  Kef.  2-1  shows  that  two  adjacent  panels 
separated  by  stringers  are  uncorreluted  by  running  waves  but  correlated  for 
low  frequency  modes  corresponding  mainly  to  those  modes  near  the  fundamental 
frequency  of  the  entire  bay  of  panels.  This  section  of  the  paper  is  directed 
toward  accurate  assessment  of  the  panel  response  to  turbulence  In  terms  of 
model  shape,  cross-spectral  density,  and  mean  square  displacement. 

h.  Theoretical  Approach 

The  correlation  function  of  the  displacement  covariance  due  to  a  random  force 
|>(xo,  yoto)P<x^.y^,i;>)  has  been  given  by  Hef.  15. 
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The  pressure  covariance  lias  been  reported  primarily  in  Hef.  5  for  broad  band 
response: 


(t  -  t*  ) 
'  o  o' 


<p(x  ,y  ,t  )p(x*  ,y'  ,t’  )> 
1  '  0  0  o'  '  o,Jo  o 


:! 

2 


*) 

=  p~  O 


^  V 


i-=l  ,.2/jr  (fx  -x’  ]  -  U  It  -t'  })”  +  (y  -  y'  )‘ 

IVy  !  ^ jTy/  I'  1  O  O'  C1  O  O1'  W0  J  o' 


A 


2 

2  if-  (16) 

e  1  K«' 


'l'ho  plate  input  response  function  g(x,y,t  |x^,y^,l'0)  was  evaluated  in  terms  of 
eigen-functions  or  orthogonal  modes  of  plate  oscillation  which  are  of  the  form 


Y  (x.y.t)  ~  o  (x.y)  exp  -  a  t  i  ia>  t 
m,n'  ’  ym,n'  1  m,n  m,n 


'l'ho  normal  mode  satisfied  the  equation: 


Bv‘lY(x,y>i)  ‘  Mac  1  =  0 


(17) 


(18) 


it  shows  explicit  division  of  the  damping  into  acoustic  and  structural  damping. 
Substituting  ( 17)  into  ( 18), 


v‘lp(x,y)  -  rm  n  0(x,y)  =-  o 


where 


(19) 


i  JIA  (a  i  u)  )“  -  (/f  sr.  +/W)(a  1  »w  ) 

p  1  „  u  v  m , n  in , iv  w  AC  *  STM  m,n  m,ir 
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(20) 
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The  impulse  response  function  for  low  damping  was  found  to  lie: 


g(x,y,t|x  ,y  ,to)  cc  2 

m.n 


I'm.VL&’Mm.'Po-yJ 


(J 


m,n 


exp  I  -  a  (1  -  t  )  sin  oi  (t  -  t  )  j 
1  |  m,n'  o'  m,n  1  o'} 


(21) 


and  the  eigenfunction  solution  to  the  homogeneous  wave  equation  of  the  plate  in 
the  ease  of  freely  supported  boundaries  is 


(ab) 


•  nurx  .  mrv 
1/2 sin  a  sin  b 


(22) 


An  attempt  has  been  made  to  evaluate  the  integral  (15)  using  two  methods  of 
integration;  the  Gaussian  quadrature  and  the  Monte  Carlo  method.  The  problem 
was  too  lengthy  to  handle  with  the  standard  interaction  technique.  After  some 
survey,  both  the  Monte  Carlo  method  and  the  Gaussian  quadrature  methods 
proved  unfeasible  because  an  estimation  of  the  maximum  value  of  the  integrand 
(over  six  dimensional  space)  is  necessary,  and  the  extreme  demands  on 
computer  time  made  it  impossible  to  proceed  with  this  typo  of  computation. 
Knginchring  judgement  suggests  that  certain  coupling  terms  may  be  eliminated  in 
the  integrand,  thus,  considerably  reducing  the  lime  of  integration.  Both  Dyer  and 
el  Baroudi  use  a  relatively  simple  input  response  for  pressure  correlation.  One  can 

deduce  in  the  ease  of  the  delta  function,  any  term  in  the  integral  which  includes 
the  product  of  two  identical  cigcnmodcs  effectively  vanishes  when  the  integration 
is  performed. 
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With  the  present  type  of  cross-correlation  function,  t he  closed  form  solution 
may  he  impossible  or  extremely  lengthy  using  numerical  computer  techniques. 
Therefore,  an  assumption  is  made  based  on  el  Burundi's  results,  that  the 
cross-modal  effect  be  eliminated.  This  reduces  equation  (15)  to 


— - — rrrr  p2  ^  ^ 
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By  changing  the  variables  of  integration  (see  Appendix  B),  integral  (25)  is 
reduced  to: 


Y(x,y,l)y(x',y',l') 


v  ^.n...<X’y)  . (X'-yl) 

■>,-  2  (^-Tlr2  1  > 

M  K»  Jlt> 


niff 

nff  «  :t 

/  m 

i'(y)  hr(x)  2 

-m/r 

Jo  l  o  P-i 

A 

-|(t  t  x) I/O 

pKj/  e 

A„K,,  e 


l(r-x)  |  /0 


1  (kuF  [(IS- v|-  V)2 


-  dx  dy  d7. 

hi  - 


NUMBER  D6-99H  -  VOL.  Ill 
REV  LTR 
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f(Z)  -  COS  Z  1  —  (-  jz|  COS  V.  1  sill  |z|  ) 


l'{y)  =  cos  j"  1  ^  (sin  y  -  y  cos  y) 


-a  x  r  u 

K(x)  -  c  ni,n  sin  a;  xi  -------  00s  uj  x 

n,,n  amn 


It  is  easier  now  to  obtain  further  information.  The  numerical  result  from 
(21)  {>ives  both  the  modal  and  broad  band  correlation. 


'I’he  normalized  correlation  function  is  defined  as 


Hy({.*I.T)  - 


Y(x.y,t)Y(x‘.y,tt') 
Y(x,y,t)“Yfx'yM’)“  l> 


and  the  normalized  cross  power  spectral  density  oi  the  panel  displacement, 


lyKj.K^.ui)  j  «• 

<«»>*/$  y  'T) l' 


i(Kj5  •  K.,») '  wt) 


il^drjdr 


The  panel  and  turbulent  spectra  are  related  by 


r(KJfK2,w)  "  -  I‘y(KrK2,«)  "  I’tKj.K^w) 


The  transfer  function  l'(K  , K.^.w)  can  now  be  easily  computed  since  both  Input 
and  output  spectra  are  known.  From  acoustic  interest  one  would  like  to 
determine  the  modal  volume  displacement  since  it  is  related  to  the  acoustic 
power. 
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Vol  •-  S  Y  /  /  0(x)  0(y>  tbctly  (28) 

'  o  o 

i 

c.  Example: 

'i’he  above  results.  Equations  (2-1)  to  (2(>)  and  (28),  are  numerically  evaluated  for 
hvo  panel  sizes,  12  x  7  x  0.  0*1  inches  and  .'!(>  x  (>-  1/2  x  0.  (H  inches,  corresponding 
to  the  panel  sizes  used  in  the  laboratory  experiment.  Results  cun  be  compared 
only  qualitatively  with  experiment,  because  the  experimental  panels  arc  clamped 
rather  then  simply  supported  as  assumed  in  the  theory. 

Let  it  be  assumed  that  U  /U  =  0.  8  and  6*  -  0.155  inch.  The  moan  square 
displacement  is  obtained  from  Equation  (25)  by  setting  r  -  o;  x  -  x ' ;  y  -  y 1 ;  l  -  tj . 
The  value  of  x.y  are  chosen  such  that  the  displacement  will  be  maximum. 


The  values  of  total  damping  for  the  12x7  x  0.  01  inch  panels,  assumed  to  be 
a  -  .5  (w  )1//<!  and  a  -  5  (o  , for  the  OG  x  G-l/2  x  0.01  inch 
panels,  has  been  measured  on  the  same  size  panel  with  the  clamped  edges  when 
excited  by  a  turbulent  boundary  layer.  One  interesting  [joint  is  that  the  modal 
damping  is  considerably  different  when  the  panel  is  excited  by  turbulent  flow 
then  when  excited  by  pure  tones. 


The  mean  square  displacement  computed  at  the  center  of  each  mode  is  shown 
in  Fig.  12  for  the  12x7x0.  01  inch  panel.  It  increases  with  the  convection 
velocity  for  the  lower  order  modes  while  for  the  highest  modes  the  mean 
square  displacement  levels  oil-  for  large  values  of  Uy.  This  result  is  quite 
consistent  with  the  previous  experimental  work  reported  by  the  author  in  Ref.  5 
which  shows  that  the  sudden  change  in  slope  is  due  to  mismatch  in  wave  number 
and  frequency  between  the  turbulent  and  panel  spectra. 


The  experimental  results  indicated  also  that  above  the  peak  displacement  the 
mean  square  amplitude  decreases  for  further  increases  in  Mach  number.  This 
pronounced  reduction  is  not  fell  by  the  theoretical  panel  because  the  decay  of 
the  eddy  lifetime  0  with  increasing  convection  velocity  U  is  not  appreciable 
since  0  is  a  weighted  average  over  a  narrow  bandwidth  (Ref.  5,  Fig.  5)  due  to 
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limited  range  of  frequency  component  of  the  correlation  measurement. 

Therefore  the  experimental  panel  favors  a  larger  mismatch  because  the  wall 
pressure  spectrum  PfK^w)  decays  faster  than  the  panel  spectrum  P^K^te).  The 
comparison  between  theory  and  experiment  for  the  12  x  ?  x  .0-10  panel  using  the 
measured  value  of  0  is  shown  in  Fig.  13.  Figure  12  shows  also  the  variation 
in  mean  square  displacement  with  total  dumping.  The  change  in  damping  con¬ 
tributed  to  an  almost  linear  change  in  the  mean  square  displacement.  A  more 
appropriate  value  ofO  .  total  damping  a|n  and  modal  bandwidth  has  been  selected 
from  experimental  obser  /alien  in  comparing  the  measured  mean  square  displacement 
with  the  correlated  one  for  the  3G  x  G.  5  x  0.  >10  inch  panel  reported  in  Figs.  14 
and  15.  Noticeable  is  the  shift  of  the  maximum  mean  square  displacement  with  U 
due  to  the  difference  in  frequency  between  the  two  panels  because  of  the  different 
edge  conditions.  The  comparison  is  satisfactory  to  validate  the  reliability  of  the  theory 
in  predicting  the  response  of  a  single  panel,  with  a  proper  choice  of  the  value  of 
the  damping  and  eddy  lifetime. 


The  normalized  theoretical  and  experimental  broadband  space  time  correlation 
for  a  3G  x  G.5  x  .010  inch  panel  is  shown  in  Fig.  16.  The  two  autocorrclograms 
taken  at  the  center  of  the  panel  are  quite  similar,  implying  that  the  number  of 
modes  contributing  to  the  broadband  response,  40  in  all,  that  is  m  =  (1  to  20) 
and  n  ^  (1 ,2)  and  the  associate  modal  damping  factors  are  close  to  those  of  the 
expert  mental  panel. 


Results  can  be  extended  to  obtain  the  cross- correlation  everywhere  along  the 
panel  by  simply  multiplying  by  the  eigenvalue  of  Kq  nation  (24). 

Due  to  the  complexity  of  the  triple  integration  (Kquation  (24)),  the  rigorous 
validity  of  the  numerical  result  may  depend  on  the  previous  assumption  that  the 
cross  modal  effect  is  eliminated;  however,  experience  obtained  from  this 
calculation  has  indicated  that  the  complete  problem  can  be  solved  including  the 
cross- modal  effect. 
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The  transfer  functions,  defined  in  Equation  (2?)  as  the  ratio  between  power 
spectral  density  of  the  turbulence  and  the  panel  in  the  wave  number  frequency 
scale,  constitute  the  most  important  result.  In  the  past,  a  similar  approach 
was  used  (Ref.  Id)  expressing  the  transfer  function  in  terms  of  frequency  and 
wave  number,  but  the  integration  proved  to  be  a  stumbling  block.  Theoretically, 
the  one-dimensional  wave  number  frequency  spectrum  of  the  exciting  pressure 
field  is  obtained  from  Equation  (G),  (computed  values  are  shown  in  Fig.  4)  and  the 
panel  response  spectrum  can  be  obtained  from  the  Fourier  transform  Equation  (2G)  of 
the  cross-correlation  function  (Equation  24).  Experimental  spectral  densities  of 
both  turbulence  and  panel  are  shown  in  Fig.  18  while  theoretical  panel  spectral 
density  lias  not  yet  been  computed.  The  effect  of  matching  and  mismatching  the 
wave  number  for  a  constant  frequency  clearly  shows  the  mechanism  of  the 
excitation  of  the  structure  by  turbulence.  The  experimental  results  have  been 
chosen  such  that  a  strong  coincidence  effect  dominates.  This  is  also  indicated 
by  the  closely  matching  behavior  of  the  two  eospoolra.  An  interesting  point 
is  that  the  turbulent  eddy  decays  much  faster  than  the  wave  on  the  panel.  There¬ 
fore,  the  matching  of  the  two  spectra  occurs  only  in  a  very  narrow  range  of 
frequency  and  wave  number. 

The  modal  volume  for  the  .'{(5  x  (i.  5  x  .  040  inch  and  :JG  x  G.  5  x  .  08  inch  panels 
has  been  computed  from  Equation  (28)  to  predict  the  acoustic  power  radiated  of 
Section  VI,  since  the  power  radiated  is  proportional  to  the  uncancclcd  volume 
velocity. 
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V.  EXPERIMENTAL  RESULTS  ON  THE  REFLECTION  AND  TRANSMISSION  OP  A 
FLEXURAL  WAVE  ON  THE  PANEL 


Space  time  correlation  (or  various  panel  displacements  consistently  indicates 
that  the  panel  response  Is  prominently  excited  by  running  flexural  waves  with 
the  highest  correlation  oeeuring  along  lines  £-  U  t  -  constant  in  the  £  ,r  plane 

I  ^ 

as  shown  in  Fig.  17  where  £-  x  -  x  .  (Ref.  5  and  21).  This  demonstrates 
that  the  average  convection  velocity  coincides  with  the  average  wave  velocity 
in  the  panel.  The  correlograins  also  indicate  that  the -surface  waves  have 
a  phase  velocity  of  ±U  The  symmetry  in  the  positive  and  negative  time  delay 
is  interpreted  as  waves,  one  moving  with  the  flow  and  another  opposite  to  the 
flow. 


The  reflection  and  transmission  problem  is  experimentally  investigated  on  a 
.'Hi  x  ti-l/2  x  0.  04  inch  panel  by  setting  up  two  test  configurations,  in  the  first 
configuration  the  panel  is  damped  at  the  downstream  side.  The  damping  mate¬ 
rial  consists  of  a  1/2- inch  thick  layer  of  sand  retained  within  4  inches  of  the 
end  by  soft  rubber  wedges  glued  across  the  panel.  The  spaoe  time  correlation 
of  the  displacement  measured  along  the  centerline  from  the  center  toward  the 
damped  side  is  shown  in  Fig.  19.  The  space  time  correlogram  for  the  panel 
with  damped  edges  is  quite  different  from  the  correlogram  for  the  bare  panel. 
The  wave  moving  with  +U  (downstream)  is  undamped  for  both  panels  while  the 
wave  moving  with  -Uc  (upstream)  the  dumped  panel  shows  that  the  negative 
correlation  is  damped.  The  interpretation  is  that  the  sand  has  partially  damped 
the  incident  wave  at  the  downstream  edge  since  no  correlation  and  convection 
results  with  the  negative  delay. 

The  second  configuration  is  made  on  the  same  si/.e  panels  divided  by  two  equally 
spaced  stringers  glued  to  the  panel,  as  reported  in  Ref.  24.  The  correlation 

1  •  • 

between  stringers, (Fig,  20)  shows  the  well  known  behavior  of  waves  moving 
with  ±UC  while  the  correlation  across  the  stringers  indicates  an  unswept  pattern 
in  the  time  scale  corresponding  to  a  standing  wave  pattern.  The  interpretation 
of  the  correlogram  is  that  running  waves  are  reflected  at  the  stringer  bound¬ 
aries  while  the  standing  waves  are  transmitted. 
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VI.  RADIATION  OF  SOUND  BY  FLEXIBLE  PANELS 

In  the  previous  sections  we  have  discussed  the  response  of  panels  excited  by 
turbulent  boundary  layers,  and  in  this  section  we  will  explore  the  radiation 
resulting  from  the  vibration  of  those  panels.  In  the  structural  response  to  tur¬ 
bulence  study,  it  was  found  that  the  structural  response  is  in  the  form  of  a  forced 
vibration  and  coincident  vibration.  If  the  structural  panel  has  a  mass  law 
behavior  (forced  vibration),  the  resixmse  will  then  be  computed  easily  since  it 
depends  on  the  weight  per  unit  area  of  the  structure;  in  the  case  of  coincidence, 
the  response  is  above  the  mass  law  and  the  problem  becomes  more  complex 
since  the  pressure  field  on  the  panel  becomes  the  controlling  mechanism  of  the 
response.  To  estimate  the  sound  power  level,  the  modal  volume  displacement 
or  mi  equivalent  radiation  efficiency  must  be  ascertained,  Theroetieal  pre¬ 
diction  of  the  radiation  efficiency  and  measurement  of  the  acoustic  damping 
have  been  reported  in  Refs.  25  through  27,  and  5.  The  radiation  property  of 
the  panel  is  related  to  the  volume  velocity  in  each  mode,  which  often  radiates 
like  independent  monopoles.  For  a  finite  panel,  the  major  source  of  radiation 
below  the  critical  frequency,  f  =  l/2ira“/KC|, arises  from  the  interaction  of  the 
bending  wave  with  the  discontinuity  ol'  the  boundary.  The  present  panel  sizes 
have  modes  that  radiate  mainly  below  the  critical  frequency  where  the  radiation  is 
is  somewhat  less  efficient  than  it  is  at  the  critical  frequency  or  above.  Below 
the  critical  frequency,  the  radiation  from  a  1/4  wavelength  of  the  mode  segment  ii 
is  cancelled  by  the  radiation  of  the  adjacent  1/4  wavelength  which  is  out  of  phase. 
This  process  of  cancellation  is  extended  across  the  four  boundaries  of  the  panel 
and  the  effective  radiation  area  lying  between  the  panel  edges  and  the  nearest 
1/4  wavelength.  If  the  flexural  wavelength  component  is  larger  than  the  acous¬ 
tic  wavelength,  as  in  the  case  of  the  17-1  and  19-1  modes,  of  the  36  x  6-1/2  x.  040 
inch  panel  (Fig.  21)  the  1/4  wavelengths  along  the  two  side  edges  are  acoustically 
uncoupled  and  combine  to  form  an  edge  mode.  Since  the  edge  mode  occurs  at 
both  ends  of  the  panel  whose  corners  are  36  inches  apart,  the  edge  mode  radia¬ 
tion  becomes  uncoupled,  because  the  distance  between  them  is  much  larger  than 
the  acoustic  wavelength.  Each  edge  mode  gives  rise  to  a  edge  radiator  which 
radiates  independently.  This  cancellation  process  can  be  carried  out  for  all  m 
and  n  modes ,  which  shows  that  the  resulting  sources  are  monopole  for  the  odd-odd 
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mode,  dipole*  for  the  even-odd  mode,  and  quadruple  for  the  even-even  mode. 
Since  power  radiated  by  the  dipole  or  quadrupole  in  much  smaller  than  the 
monopole,  one  can  neglect  the  power  in  the  even-odd  and  even- even  mode.  In 
this  analysis,  the  modal  volume  displacement  has  been  measured  for  the  odd- 
odd  and  for  the  even-odd  modes  associated  with  the  acoustic  power  radiated  in. 
a  reverberation  room.  Both  measurements  are  useful  in  making  an  accurate 
assessment  of  the  radiated  power  to  be  used  with  the  response  theory  of  the 
panel  displacement  of  the  previous  section  to  predict  the  acoustic  power  radiated. 
The  modal  acoustic  power  radiated  in  a  reverberant  field  in  terms  of  modal 
velocity  is  given  by 


NoTPeK  2— rr  2P  +  P 
PVVL  =  - —Y1  — — l 


m .  n 


4tt 


Kwl  l*..wl Ilxdy 


•2* 
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The  value  of  N  depends  on  the  radiation  mode  which  can  be  classified  as  an  edge 
radiation  or  piston  radiation.  Here,  piston  radiation,  for  which  N  -  4,  was 
assumed  for  the  purpose  of  eaieulation.  The  value  of  N  also  depends  on  the  distance 
between  sources  in  relationship  to  the  acoustic  wavelengths  as  well  as  on  the 
panel  Iwundary  (Refs,  25  and.  27).  The  ratio  (21^,  t  P  )/P  is  the  contribution 
due  to  stringers  mounted  on  the  panel,  following  Ref.  25.  Comparison  of  the 
radiated  power  is  shown  in  Figs.  22,  23  ,  and  24  by  direct  measurement  and  by 
Equation  (50)  using  a  measured  value  of  the  mean  square  displacement.  Thu 
comparison  is  in  satisfactory  agreement  for  a  simple  panel  for  thicknesses 
.  080- inch  and  .  040- inch  as  well  as  for  panels  with  stringer  combinations  using 
the  idealized  ratio  suggested  by  Ref.  25.  It  is  believed  therefore  that  with  the 
proper  choise  of  damping,  the  total  acoustic  power  is  obtainable  directly  from 
the  structural  response  of  the  panel  given  by  Equation  (24). 

1 

The  eigenvalues  for  (he  lateral  and  longitudinal  modes  of  the  3 G- inch  by  G.  5-inch 
panel  with  rigid  boundaries  are  given  in  Table  I,  following  the  same  method  used 
in  Refs.  5  and  29, 


*  The  use  of  this  equation  has  been  suggested  to  the  author  by  Dr.  R.  II.  Lyon. 

. . - — _ _ _ _ _ 
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Another  interesting  method  for  predicting  the  radiated  power  has  been  recently 
reported  by  Ref.  23  using  the  modal  density  approach  averaging  over  a  band¬ 
width  for  vibration  response  and  radiation.  This  approach  shows  good  agree¬ 
ment  with  two  different  experiments,  one  of  which  is  that  made  by  the  author. 

The  most  relevant  phenomena  relevant  to  both  power  radiation  and  the  response  of 
the  panel  is  the  effect  of  coincidence,  characterized  by  the  peak  in  amplitude  dis¬ 
placement,  aswellas  in  the  acoustic  power  radiated.  Above  the  peak  displacement, 
the  effect  of  rolling  off  in  amplitude,  appears  to  be  associated  with  progressive 

decrease  in  wave  length  which  contribute  to  the  change  in  the  power  law  from 
5  2  3 

U  to  IT’  discussed  in  the  previous  experiment  Refs.  4  and  5.  This  change  can 
be  easily  visualized  from  the  cross-power  spectrum  of  both  turbulence  and  panel 
(Fig,  18).  The  mismatch  both  above  and  below  coincidence  in  frequency  and 
wave  number  is  the  cause  of  the  changing  power  law  and  can  be  attributed  to  a  rapid 
decay  of  the  eddy  in  the  turbulent  boundary  layer  in  relationship  to  the  panel  mode. 

VII.  RADIATION  OF  SOUND  FROM  HONEYCOMB  PANELS 
The  honeycomb  structure  has  recently  received  wide  application  by  airframe 
manufacturers  because  of  its  Integrity  and  light  weight  characteristics.  There¬ 
fore,  it  became  useful  to  determine  the  noise  radiated  from  such  a  structure 
under  boundary  layer  excitation.  Measurements  were  conducted  on  a  size 
7  x  12  inch  panel  with  3  types  of  honeycomb:  two  having  total  thickness  of 
0.  25  inch.  One  panel's  surfaces  was  made  of  0.  01  inch  aluminum;  the  other 
with  0.  01  inch  titanium.  The  third  honeycomb  panel  was  0. 1  inch  thick  with  the 
surfaces  made  of  an  0.  01-inch  aluminum  sheet.  All  three  panels  have  an  aluminum 
core  glued  to  the  top  and  bottom  panels.  Since  all  three  panels  have  a  weight 
comparable  to  the  common  aircraft  panel,  a  0.  040- inch  thick  aluminum  panel  is 
chosen  for  comparison  of  the  acoustic  power  level  rutilated  using  a  rudiating 
area  of  7  x  12  inches.  The  acoustic  power  is  plotted  in  Fig.  25  in  l/3- octave 
bands.  The  level  of  the  0.  040- inch  panel  dominates  over  the  honeycomb  panels 
with  lower  frequency  characteristics.  The  second  highest  level  is  the  aluminum 
honeycomb  0.  1-inch  thick  which  shows  a  shift  in  the  fundamental  mode  from 
300  to  000  cps  with  respect  to  the  0.  040  panel.  A  considerably  lower  level 
results  from  the  0.  25- inch  aluminum  and  titanium  panels  which  show  a  com¬ 
parable  level  between  them  except  for  the  3  db  differences  at  the  peak.  For  the 
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aluminum  panel  the  peak  occurs  at  J000  cps  and  for  the  titanium  at  the  1250  cps. 

A  slight  shift  in  the  maximum  response  is  expected  due  to  the  slight  change  in 
weight  between  the  aluminum  panel  {0.  33  pounds)  and  the  titanium  panel 
(0.  26  pounds). 

I 

The  structural  response  of  the  honeycomb  panel  has  not  yet  been  determined  but 
the  preliminary  results  indicate  that  the  response  predominates  at  a  higher  fre¬ 
quency  than  with  the  simple  aluminum  panel.  The  highest  frequency  sound 
radiation  of  the  honeycomb  can  be  very  easily  attenuated  by  the  usual  absorbent 
acoustic  material  whereas  in  the  lower  frequency  region,  the  simple  aluminum 
panel  is  still  a  major  problem  for  the  transmission  of  sound.  One  can  conclude 
from  this  experiment  that  the  transmission  loss  through  the  honeycomb  structure 
is  a  function  of  stiffness  which  is  related  to  thickness.  Figure  26  also  shows 
the  comparison  of  the  change  in  power  level  with  Mach  number.  At  high  speed, 
the  honeycomb  radiates  more  efficiently,  however,  the  level  is  still  lower  than 
the  usual  skin  panel.  In  die  case  of  aircraft  sidewall  structure,  a  liner  is  used 
which  satisfies  a  double  purpose:  thermal  and  acoustic  insulation.  This  gives 
considerable  sound  attenuation  in  the  higher  frequencies  and,  therefore,  improve¬ 
ment  in  level  between  the  honeycomb  panel  and  the  normal  sidewall  panel 
structure. 

VIII.  CONTRIBUTION  OF  T11E  PANEL  BOUNDARY  TO  SOUND  TOWER 
RADIATION 

Experimental  results  have  indicated  that,  for  frequencies  below  the  critical 
frequency,  the  power  radiated  by  the  punel  changes  with  the  edge  conditions 
of  the  boundaries.  Analyses  made  in  Refs.  27  and  33  huve  indicated  some 
sensitivity  of  the  rudiation  efficiency  to  the  edge  conditions.  The  results  com¬ 
puted  in  the  previous  chapter  confirmed  that  the  noise  comes  from  the  unean- 
celled  edge  half  because  of  the  interaction  with  the  boundary.  A  series  of 
experiments  have  been  conducted  with  different  edge  boundaries  to  determine 
if  there  is  any  possibility  of  reducing  the  acoustic  power  radiated.  An  interest¬ 
ing  result  obtained  demonstrates  the  possibility  of  minimizing  the  acoustic 
power  radiated  by  modifying  the  panel  boundaries.  The  best  results  show  an 
uverage  of  10  db  reduction  in  the  power  level  by  using  15  percent  additional 
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panel  weight  to  modify  the  boundaries.  Typical  data  using  three  boundary  con¬ 
figurations  on  a  30  x  0-1/2  x  0.040  inch  panel  are  shown  in  Figs.  27,  28,  and  29. 

The  first  configuration  is  made  up  of  three  layers  of  different  lengths  of  0.001 
inch  thick  steel  membranes  glued  to  the  upstream  and  downstream  edges  of  the 
panel.  Ettectively,  the  membranes  increased  the  panel  stillness  gradually 

t 

through  the  last  wavelength  toward  the  boundaries. 

The  second  configuration  consists  of  rubber  wedges  at  both  ends  of  the  panel. Fig.  30 
The  rubber  wedges  are  lighter  than  the  membranes  with  an  optimum  reduction 
in  acoustic  power  above  250  cps.  Figures  27  and  28  show  the  results  in  terms 
of  flow  Mach  number  compared  with  a  bare  panel  and  a  panel  mounted  on  the  two 
equally  spaced  stringers  and  Ftg.  29  shows  the  results  in  l/3-octave  band  level. 

The  structural  response  of  the  first  two  configurations  has  not  yet  been  made. 
However,  it  is  expected  that  an  increase  in  the  frequency  and  decrease  in  the 
volume  displacement  will  occur.  This  last  decrease  will  contribute  a  lower 
radiation  level.  Speculation  can  be  made  that  a  reduction  in  level  may  have 
some  contribution  from  the  mismatch  in  the  frequency  and  wave  number  scale 
from  waves  propagating  toward  the  boundary  of  the  panel  due  to  low  damping 
and  from  the  wave  length  resulting  from  the  corner  interacting  with  the  phase 
and  amplitude  of  the  local  turbulent  mode.  The  reduction  in  power  level  is 
significant  enough  to  encourage  further  investigation. 

The  last  configuration  consists  of  damping  the  flexural  wave  on  the  panel  propa¬ 
gating  with  the  speed  ±U  .  The  downstream  side  of  the  panel  only  is  damped  by 
sand  to  interpret  more  significantly  the  reflected  part.  This  is  the  space  time 
correlation  case  discussed  in  Section  V  and  in  Figs.  17  and  19  with  and  without 

damping.  The  waves  that  are  not  reflected  from  the  downstream  edge  of  the 

» 

boundary  contribute  significantly  to  a  reduction  in  acoustic  power.  This  method, 
though  impractical,  shows  the  contribution  of  the  running  wave  to  the  acoustic 
power  radiated  when  it  is  reflected  by  the  boiunlary. 
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IX.  EFFECT  OF  PRESSURIZATION  ON  THE  RESPONSE  OF  THE  PANEL 
EXCITED  BY  A  TURBULENT  BOUNDARY  LAYER 

An  aircraft  flying  at  high  altitude  has  a  static  pressure  differential  between  the 

cabin  pressure  and  the  local  atmosphere.  This  produces  a  static  deflection  on  the 

fuselage  panels  which  is  related  to  the  geometry  and  framework  of  the  structure 

and  loading.  If  the  deflection  is  not  small  compared  to  the  thickness  of  the 

panel,  the  stretching  of  the  middle  surface  of  the  plate  must  be  taken  into 

account  in  calculating  the  modal  frequency  of  vibration.  Duo  to  the  resistance 

of  the  plate  to  stretching,  the  rigidity  and  frequency  will  increase  with  the 

pressure  differential.  Measurements  made  in  the  boundary  layer  facility  have 

shown  that  the  frequency  of  vibration  and  the  radiation  level  shift  to  higher  fro- 

t 

queney  with  the  increase  in  pressure  differential  across  the  panel.  The  correct 
magnitude  of  the  displacement  level  and  acoustic  power  level  has  not  boon 
validated"  because  the  curvature  of  the  flexible  panel  joining  the  rigid  one  at  each 
of  the  four  boundaries  will  cause  the  flow  to  trip  and  generate  additional  noise. 

The  sudden  change  in  slope  as  the  flow  approached  the  larger  static  deflection  of 
the  boundaries  which  has  superimposed  on  dynamic  vibration  is  a  source  of 
additional  noise. 

Tests  on  the  pressurization  effect  were  possible  on  the  720  airplane  both  for  the 
acceleration  level  on  a  sidewall  panel  and  sound  pressure  level  taken  a  few  inches 
away  from  the  side  panel.  Results  shown  in  Fig.  31  reasonably  indicate  a  shift 
due  to  change  In  pressurization  from  3  to  8  psi  differential  at  an  altitude  of 
25,000  feet  and  Mach  0.  87,  The  change  in  amplitude  of  the  acceleration  level 
with  the  pressurization  cannot  be  regarded  as  a  general  increase  since  one  has  to 
obtain  the  space  average  of  the  mode  rather  than  the  rms  value  at  a  point.  The 
frequency  shift  due  to  pressurization  results  of  the  order  of  the  pressure  ratio 
between  the  cubing  and  the  pressure  attitude.  A  freise  prediction  is  given  in 
Ref.  80  where  agreement  with  the  measurements  is  excellent. 

X.  CONCLUSIONS 

Measurement  and  analysis  of  the  statistical  properties  of  the  surface  pressure 
fluctuation,  the  response  of  panel,  aiul  the  resultant  radiation  field  were  carefully 
Investigated  and  its  most  important  feature,  the  response  and  acoustic  radiation, 
can  be  reasonably  predicted  given  the  initial  flow  condition  and  panel  damping. 
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The  wave  number  and  frequency  spectrum  of  the  pressure  fluctuation  show  that 
the  main  energy  is  eoneentraf  around  the  convection  velocity  with  the  highest 
level  at  the  lowest  frequency  at  ive  number,  Indieating  that  eaeh  peak  is  eon- 
vectcd  at  different  velocities.  For  a  given  frequency  a  critical  wave  number 
exists,  below  which  the  wall  pressure  is  part  of  an  oscillatory  field  or  sound 
field,  and  the  eddy  convection  velocity  becomes  supersonic. 

Comparison  has  been  made  of  the  cross  power  spectrum  of  the  turbulence  with 
the  cross  power  spectra  of  the  panel.  The  superposition  of  wave  number  and 
frequency  spectra  indicate  the  matching  or  mismatching  of  spectra,  the  most 
significant  excitation  of  the  panel  by  turbulence. 

The  acoustic  power  radiated  comes  from  the  uncancelled  volume  velocity  at 
the  panel  boundary  in  agreement  with  Ref.  25.  By  modifying  the  boundary 
condition,  considerable  reduction  in  acoustic  power  radiated  can  bo  obtained. 
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APPENDIX  A 


Equation  (1)  shows  that  the  decay  of  the  broadband  space-time  correlation  is 
given  by  exp  (-|£|/UC0)  rather  than  the  more  common  form  shown  in  previous 
reference  exp  (-|t|/0).  It  is  rather  difficult  to  see  the  implication  of  the  choice 
of  the  exponential  without  performing  the  transform  with  respect  toT,  as  well  as 
by  comparing  the  cross  correlation  using  both  exponential  forms.  Therefore, 
rewriting  Equation  (1)  using  the  form  exp  (-|r!/0)  we  obtained: 

-M  /Q  -|(«-U  ■ T)2+H2|/2a12  r>2+  22j/2cr.2 

A.  e  1  '  +A..  e  1  c  '  1  (!') 


Taking  the  Fourier  transform  respect  to  r  Equation  (P)  becomes 
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The  power  spectrum  on  (2')  is  of  the  same  type  as  the  one  given  by  Equation  (13) 
using  exp  -|£|/UC0  form.  The  differences,  however,  are  in  the  phase  change, 
by  the  factor  a-/Ec^0  which  is  constant.  The  comparison  of  the  cross  correla¬ 
tion  is  given  by  Equations  (1)  and  (J')  as  well  as  shown  in  Fig.  32. 

Since  excellent  agreement  indicates  the  functional  form  of  the  wall  pressure 
correlation  fits  the  experimental  cross  correlation  distribution,  either  one  ol 
two  lorms  of  exixmential  may  be  used  to  describe  the  decay  ot  the  wall 
pressure  correlation. 
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APPENDIX  LI 


DERIVATION  OF  THE  CROSS-CORUE  RATION  -  INTEGRAL  21 


In  order  to  simplify  the  computation  of  Equation  (24),  an  attempt  is  made  to 
change  the  variable  of  integration,  thus  reducing  the  number  of  operations  for 
the  Integral.  Rewriting  Equation  (24) 
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The  mode  shape  which  satisfies  the  eigen  function  equation  subject  to  the  panel 
boundary  condition  for  simply  supported  edges  is 
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Considering  tins  .space  integrals  in  (ill)  and  using  the  product  in  (2D) 
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Similarly  for  longitudinal  coordinates 
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Using  the  definition  of  4B 

jf  f(x)  dx  -  f  (g(t))  |jy(t)!dt 
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Let  the  time  Integral  be  simplified  by  changing  the  limit  of  integration. 
From  (11B),  Ij, 
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